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We study resonant x-ray scattering (RXS) at Np M4,5 edges in the triple-k multipole ordering 
phase in Np02, on the basis of a localized electron model. We derive an expression for RXS 
amplitudes to characterize the spectra under the assumption that a rotational invariance is preserved 
in the intermediate state of scattering process. This assumption is justified by the fact that energies 
of the crystal electric field and the intersite interaction is smaller than the energy of multiplet 
structures. This expression is found useful to calculate energy profiles with taking account of the 
intra-Coulomb and spin-orbit interactions. Assuming the Fg-quartet ground state, we construct the 
triple-k ground state, and analyze the RXS spectra. The energy profiles are calculated in good 
agreement with the experiment, providing a sound basis to previous phenomenological analyses. 



PACS numbers: 78.70.Ck, 75.25.-l-z, 75.10.-b, 78.20.Bh 



I. INTRODUCTION 

Resonant x-ray scattering (RXS) technique has at- 
tracted much attention to study spin and orbital prop- 
erties of 3d transition-metal compounds. RXS at the K- 
edge is described by a second-order optical process that 
an incident x-ray excites a Is core electron to unoccupied 
4p states and then the Ap electron is recombined with the 
core hole with emitting x-ray in the dipole process {El). 
It became widely known after the observation of intensi- 
ties on orbital-ordering superlattice spots at Mn K-edge 
in LaMnOai At the earlier stage, the spectra were in- 
terpreted as a direct observation of orbital ordering!^ 
However, subsequent theoretical studies based on band 
structure calculations revealed that the spectra are a di- 
rect reflection of lattice distortion(2i^ since the Ap state 
in the intermediate state is influenced not by the orbital 
ordering of 3d electrons but by lattice distortion through 
the hybridization with the 2p state at neighboring oxygen 
sites. 

Different from transition-metal compounds, M4.5 edges 
are available for forbidden reflection Bragg spots in ac- 
tinide compounds.Si^ The RXS spectra are more di- 
rectly reflecting multipole orderings of 5/ states, since 
the El process involves a transition from the 3d-core to 
5/ states. Each actinide atom usually carries local multi- 
pole moments, which can order at low temperatures due 
to intersite interactions such as exchange interactions. 
For such localized electron systems, RXS amplitudes are 
given by summing up contributions at each site. The 
crystal electric field (CEF) and the intersite interaction 
can be safely neglected in the intermediate state, because 
they are much smaller than the intra-atomic Coulomb 
interaction. Therefore, it may be reasonable to assume 
that the intermediate state preserves the rotational in- 
variance. Under the assumption, we derive an expression 
for the RXS amplitude in the El process to character- 



ize the spectra. Although the expression is essentially 
the same as the formula by Hannon et al.^ the present 
form is useful to calculate energy profiles with taking full 
account of multiplet structures. Using this expression to- 
gether with a microscopic model, we calculate the RXS 
spectra in the triple-k multipole ordering phase in Np02. 

Np02 undergoes a second-order phase transition below 
To = 25.5 K.^°'^^ Since Np ions are Kramers ions in the 
(5/)'^ configuration, a magnetic ground state is naturally 
expected. However, neither Mossbauer spectroscopj*iSii^ 
nor neutron diffraction experimentsi^*i^ could detect any 
evidence of the sizable magnetic moment. Actually, the 
former experiment gave an estimate of the upper limit of 
the magnitude of the magnetic moment ~ 0.01/iB, which 
was too small to explain the effective paramagnetic mo- 
ment ~ 2.95/xb^ Another complication is that a muon 
spin relaxation (/j,SR) experiment has suggested the low- 
temperature phase of breaking time- reversal symmetry^ 

A natural way to reconcile with the above observations 
is to introduce the higher-rank multipole ordering rather 
than the dipole moment. Actually, Santini and Amoretti 
proposed a octupole ordering of r2 (xyz) symmetry^Siifi 
However, this phase can be ruled out because it gives 
rise to no RXS intensity. Recently, Paixao et al. have 
reported that a longitudinal triple-k octupole ordering 
accounts well for their RXS experiment ^24 The reason 
for anticipating triple-k orderings is that it excludes a 
crystal distortion or a shift of oxygen positions, which 
is consistent with the experiment. Experimental data 
obtained from the ^^O NMR spectrum, which indicate 
the existence of two inequivalent oxygen sites, support 
the occurrence of the triple-k octupole ordering phased 
Some theoretical works also have lent support to realiza- 
tion of this type of the phasei^SiS 

Assuming the Tg-quartet ground state, we explicitly 
construct a triplc-k octupole ordering state. This state 
is found to simultaneously carry a finite quadrupole mo- 
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ment, which generates the RXS intensity. Since the RXS 
amphtudes are characterized by three terms, the scalar, 
dipole, and quadrupole ones, it is not necessary to as- 
sume the existence of the hexadecapole moment instead 
of the quadrupole moment. We calculate the energy pro- 
files with taking full account of multiplet structures in 
the intermediate state. We obtain spikes-like curves at 
Np M4 edges for smaller values of the core-level width F 
as a reflection of multiplet structures. They are found to 
merge into a single peak with F > 1 eV. The energy pro- 
file with F ~ 2 eV seems to agree with the experiment. 
The azimuthal-angle dependence of the RXS spectra is 
obtained in agreement with the previous analysis i22iSi 
The present analysis provides a sound basis to the previ- 
ous phenomenological analysis. 

The present paper is organized as follows. In Sec. ^ 
we present an expression for the RXS amplitude, which 
is useful to calculate the energy profiles. In Sec. IIIII we 
analyze the RXS spectra in the triple-k octupole ordering 
of Np02 on the basis of a localized electron model. Sec- 
tion ^3 is devoted to concluding remarks. In Appendix, 
we derive the general expression of RXS characterizing 
energy profiles. 



II. THEORETICAL FRAMEWORK OF RXS 



scattered photon 




FIG. 1: Geometry of the RXS experiment. Photon with 
polarization a or tt is scattered into the state of polarization 
a' or tt' at the Bragg angle 6. The azimuthal angle tp describes 
the rotation of the sample around the scattering vector G. 



where the dipole operators x^.j^s are defined as xi^ = 
Xj, X2j = Vj, and x^,j — Zj in the coordinate frame 
fixed to the crystal axes with the origin located at the 
center of site j. The |V'o) represents the ground state with 
energy i?o, while |A) represents the intermediate state 
with energy Ea ■ The F describes the life-time broadening 
width of the core hole. 



A. Second-order optical process 



In the resonant process, an incident photon with en- 
ergy ftw, wavevector k, and polarization vector e excites 
a core electron to an empty valence shell of the interme- 
diate state, then the excited electron falls into the core 
state emitting a photon having the same energy, wavevec- 
tor k', and polarization vector e'. For example, at M4 5 
edges in actinide compounds, a 3d core electron is pro- 
moted to partially filled 5/ states at each site by the El 
transition. The definition of a geometrical arrangement 
adopted here is found in Fig. ^ The RXS amplitude is 
assumed as a sum of contributions from individual ions. 
Since the dipole matrix element involves well-localized 
wavefunction of core states, the assumption seems quite 
reasonable. Accordingly, the RXS intensity observed in 
the experiment may be expressed for the scattering vec- 
tor G {= k' - k) as 



/(e', e, G, oj) oc 



'M,(e',e,a;) 



(2.1) 



where Mj{e' , e, oj) represents the RXS amplitude at site j 
with N being the number of sites. For the El transition, 
it is expressed a o^'^^i^^'^'' 



Mj{e', e, ui) 



(V'o|2:a,j|A)(Aka'jlV'o) 



fii^ - (Sa - E^) 



■iV 



(2.2) 



B. Energy profiles 



In localized models, the ground state and the interme- 
diate state at each site are well specified by the eigenfunc- 
tions of the angular momentum operator, |J, to). The 
CEF and the intersite interaction usually lift the degen- 
eracy in the ground state. Thus the ground state at site 
j may be expressed as IV'o):; = '^)- On the 

other hand, in the intermediate state, we can neglect 
the CEF and the intersite interaction in a good approx- 
imation, since their energies are much smaller than the 
intra-atomic Coulomb interaction and the spin-orbit in- 
teraction (SOI) which give rise to the multiplet structure. 
Thus the intermediate state preserves the rotational sym- 
metry. Under the assumption, as derived in Appendix, 
we obtain a general expression of the scattering ampli- 
tude at site y. 



Mj(e', 6,0-1) = 00(0-1)6' • e 

- zai(w)(e' X e) • (t/iqIJI^o) 

^P,(e',e)(^Aok.|V'o), (2.3) 
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where 



and 



Zl = 


Qx^ — 


2 


Jy), 


(2.4a) 


Zl = 






- J{J + 1)), 


(2.4b) 


z-i = 


Qyz 


2 ^ 


~ J zJy^ f 


(2.4c) 


Z4 = 


Q zx 






(2.4d) 


25 = 


Qxy 






(2.4e) 




e) = 


V3 , 


^y^y)j 


(2.5a) 




e) = 


- 2(24^--' 


'-x^x ~ ^y^y)l 


(2.5b) 




e) = 


V3, , 

= — + 




(2.5c) 


Piie' 


e) = 


V3 , 
= — 




(2.5d) 




e) = 


V3 , 




(2.5e) 



Here we have suppressed the dependence on j in the right 
hand side of Eq. (|2.3|l . The energy profiles are given by 
only three functions, q;o(w), ai{uj), and Q!2(w), whose 
expressions are explicitly given in Appendix. 

Several facts are immediately deduced from Eq. H2.3|l . 
First, since the scalar, dipole, and quadrupole terms ex- 
haust the amplitude, the octupole ordering alone does 
not give rise to the RXS amplitude. Second, the choice 
of the CEF parameters in the ground state does not affect 
the shape of energy profiles ao(w), ai{w) and a2('^), al- 
though it affects the expectation values of dipole and/or 
quadrupole operators. Third, aQ{ijj) has no contribution 
to the forbidden Bragg spots in the antiferro-type struc- 
ture. In order to calculate the energy profiles, however, 
we need to know explicitly wavefunctions of the interme- 
diate state, which are discussed in the next section. 



C. Absorption coefficient 

Within the El transition, the absorption coefficient is 
given by 



j a A 



r/TT 



(2.6) 

where |A) with energy E^ represents the final state, 
which is equivalent to the intermediate state of RXS. A 
comparison of Eq. (|2.t)|l with Eq. H2.2|l leads to 

A(w) cx -ImQ!o(w), (2.7) 

where Im X denotes the imaginary part of X. 




: Np 
:0 



FIG. 2: Crystal structure of Np02. Gray circles denote Np 
ions and open circles represent O ions. 



TABLE I: Slater integrals and spin-orbit interaction param- 
eters in the (3d)^''(5/)'' configuration within the HF approx- 
imation (in units of eV)i2^ 



i^'=(3d,3d) 


F*(3d,5/) 


F\^f, 5/) 


G'(3d,5/) 


180.1 






19.61 




F"^ 92.04 




f2 


9.909 




59.28 




F^ 


6.491 








F' 


4.769 




Cad = 76.254 


Cs/ = 0.298 









III. RXS SPECTRA FROM Np02 
A. Quartet ground state 

Np02 has the CaF2 type structure (F^-^j^) with a 
lattice constant a = 5.431 A at room temperature, as 
schematically shown in Fig. [SJiS Np ions are tetravalent 
in Np02, as confirmed by the isomer shift in Mossbauer 
spectra^^ and by the neutron diffraction experiment .S^ 
In a localized description, each Np ion is in the (5/)'^- 
configuration. The Hamiltonian of Np ions consists of 
the intra-atomic Coulomb interaction between 5/ elec- 
trons in addition to the SOI of 5/ electrons. The Slater 
integrals for the Coulomb interaction and the SOI pa- 
rameters are evaluated within the Hatree-Fock approx- 
imation (HFA))2S and are listed in Table ^ Because 
the isotropic parts of the Coulomb interaction F^'s are 
known to be well screened in solids compared to those of 
the anisotropic parts, the former quantities are multiplied 
by a factor 0.25 while the latter's are by 0.8. Within the 
HFA, the ground state has the ten-fold degeneracy corre- 
sponding to J = 9/2 multiplet. The choice of the multi- 
plying factors does not alter this conclusioui^ Note that 
these states of J — 9/2 are slightly deviated from those of 
the perfect Russell-Saunders (RS) couphng scheme with 
L — 6 and S — 3/2 due to the presence of the strong SOI. 
For instance, and take values 39.752 and 3.237 re- 
spectively, compared to the RS values 42 and 3.75. 

In crystal, the ten-fold degeneracy is lifted by the CEF. 



Under the cubic symmetry, the CEF Hamiltonian Hqef 
may be expressed as 



CEF 



Bi{0° + 50i) + B6{0°~210t 



(3.1) 



)^.'s represent Stevens operator equivalence. 
Thereby the degenerate levels are split into one dou- 
blet Tq and two quartets rg^-* and t'^\ The level 
scheme has been analyzed by the inelastic neutron scat- 
tering, which yields an estimate of CEF parameters as 

B4 = -3.03 X 10-2 meV and Bq = 2.36 x lO""* meV.^ 

(2) 

The lowest levels are given by the Fg , which is separated 

about 55 meV from another quartet Fg"'^'' . Diagonalizing 
Eq. H3.1|l . we obtain the bases of the lowest quartet as 



T) 
i) 
T) 

i) 



Cl 



Cl 



C4 



= C4 



C2 



C2 



C3 



C3 



(3.2) 
(3.3) 
(3.4) 
(3.5) 



with Cl = 0.2757, C2 = -0.4483, C3 = 0.8503, C4 = 
—0.9751 and C5 = 0.2216. State |m) denotes the eigen- 
state with Jz = m. Symbols t (= ±) and a (=t,i) are 
introduced to represent the state |t, cr), which distinguish 
non-Kramers' and Kramers' pairs, respectively. 




Longitudinal 




Transverse 1 



Transverse 2 



FIG. 3; Triple-k antiferro arrangements, (a) a longitudi- 
nal pattern and (b) two transverse patterns (1 and 2). Ar- 
rows illustrate the direction of the multipole moment and 
number 1, 2, 3 and 4 specify the sublattices. Oxygen 
ions are omitted. Arrows indicate vector {{Jx), {Jy), {Jz)), 
{{Qy.)AQz.),{Q^v)), and {{Ti),{T^),{T,^)), corresponding 
to the dipole, quadrupole, and octupole orderings, respec- 
tively. 



1. octupole ordering 



B. Triple-k structure 



(2) 

The four-fold degeneracy in the ground Fg ' quartet 
may be lifted by the intersite interaction, giving rise 
to induced multipole moments. Actually, several ex- 
periments tell us that the time-reversal symmetry is 
broken with nearly zero dipole moment in the ordered 
phase below Tq = 25.5 KiiLSi These observations lead 
Santini and Amoretti to propose the antiferro ordering 
of T,y,-type {T,yz = ^J^JyJ,)^ Here the over- 
line on operators means symmetrization, for example, 
JxJl = J'xJy + JyJxJy + JyJx'^^- Unfortunately, this 
phase would not give rise to the RXS intensities observed 
in the experiments. 

An important observation is that no external distortion 
from cubic structure exists in the ordered phase, that is, 
the unit cell remains cubic below Tq. This leads us to 
consider the triple-k ordering, since it allows the crystal 
to keep the cubic symmetry. As schematically shown in 
Fig. O (c), the triple-k structure is defined by all three 
members of the star of k = (001) simultaneously present 
on each site of the lattice; there are four sublattices 1, 
2, 3 and 4 at (0, 0, 0), (i, i, 0), (0, i, i) and (^,0, i), re- 
spectively. 



We start by the octupole ordering of F5„-type proposed 
by Paixao et al^ The corresponding octupole operators 
are defined by 



y 



6 

15 
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Defining operators. 



Jx{Jy — Jz), 



^(ri* + T^^ + Tf) for p= 111, 

^ {Ti - TP - TP) ior p = 111, 

^{^Ti + TyP-TP) for p= 111 

^(-Tf-Tf-KTf) for p 



= 111, 
= 111, 



(3.6a) 
(3.6b) 
(3.6c) 



(3.7) 



we assign them to each sublattice. Each Tp operator has 
eigenvalues ±ti {ti — —6.102) and doubly degenerated 
0. The eigenstates of eigenvalues ±ti are connected to 
each other by the time-reversal operation, and so are two 
degenerate states of eigenvalue 0. For example, the eigen- 
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state of eigenvalue —ti for Tm is explicitly written as 



= ^e'(«-^-t)[|+T) + e"'^|+i)] 



|-T)+e'^|-i) 



(3.8) 



with 0111 being an angle between vector (1,1,1) and 
the z axis, that is, cos 61111 = \/l/3, sin^m = ^J^J?,. 
This state is different from a ground state assumed by 
Lovesey et al, who considered the state deviating from 
Fg quarteti^ Using the eigenstates as bases, Tp is repre- 
sented as 



-h 
+ti 
000 
V 000 



(3.9) 



The ground state is given by assigning either of eigen- 
states of ztti to each sublattice; which eigenstate is rele- 
vant depends on the sign of acting mean field. As shown 
in Fig.Ol one longitudinal order and two transverse orders 
are possible in the triple-k ordering; for the longitudinal 
one, eigenstates of Tm, Tj^^^^, T^j^ and T^^ij are assigned 
to sublattices 1, 3, 4 and 2, respectively; for two trans- 
verse orders, eigenstates of of Tm, Tyii' ^iTi' ^^'^ ^iiT' 
are assigned to sublattices 1, 2, 3 and 4, and to 1, 4, 2 
and 3, respectively. 

Introducing the quadrupole operators, 



1 

V3 ^"^y 



{Qyz ~^ Q zx ^~ Qxy) 
Q zx Qxy^ 



for p ~ 111, 
for p — 111, 



Qzx - Qxy) ior p = 111, 

Qzx + Qxy) for P = 111, 



(3.10) 



we can construct the quadrupole ordering state by assign- 
ing them to each sublattice in the same way as for oc- 
tupole orderings. Since Qp's and Tp's are simultaneously 
diagonalized because of commuting with each other, Qp 
could be represented as 



Qp — 



/ -qi 

-qi 

+qi 

V +qi 



(3.11) 



with qi = -8.273. 

Let the octupole ordering be primarily realized. Then, 
each Np ion is in the eigenstate of the eigenvalue —ti (or 
ti). Since the state is also the eigenstate of the eigenvalue 
—qi, the quadrupole ordering is simultaneously induced. 
On the other hand, if the quadrupole order is primary, 
each Np ion is in the eigenstate of the eigenvalue —qi 
or qi. For the case of eigenvalue -~qi, two eigensates 
are to be degenerate and give eigenvalues —ti and ti to 
the octupole moment Tp, and thereby the net octupole 
moment becomes zero. For the case of qi, two eigenstates 
are also to be degenerate and give the eigenvalue to Tp. 
In either case, the quadrupole order carries no octupole 
order. 



2. dipole ordering 

Although the dipole ordering is ruled out in Np02, it 
may be interesting to discuss here what happens in the 
dipole ordering. Introducing the dipole operators. 




for p = 111, 

for p ~ 111, 

for p — 111, 

for p = 111, 



(3.12) 



we can construct the dipole ordering state by assigning 
them to each sublattice in the same way as in the oc- 
tupole ordering. Note that Jp and Qp are simultaneously 
diagonalized, because both operators commute with each 
other. Within the bases of simultaneous eigenstates of Jp 
and Qp, the relevant operators are represented as 





















+ji 








( 
























+32 




-91 
















-qi 








( 








+qi 
















+qi 


/ 





~h 


o\ 






-h 










V 



















0/ 





(3.13) 



(3.14) 



(3.15) 



where ji — 3.27, j2 — 0.18 with parameters given in 
Np02. The magnetic moment is evaluated on either of 
eigenstates of ±ji: {Lp + 2Sp) — 2.48 {Lp and Sp are 
defined as in the same way as Jp). 

In the dipole ordering, the ground state is given by as- 
signing one of the eigenstates of Jp's to each sublattice. 
Since ji is much larger than j2, the ground state is likely 
to be either of eigenstates of ±ji. It is obvious from 
Eqs. (|3.14|l and (|3.15|l that the dipole ordering induces 
the quadrupole moment but no octupole moment. Note 
that, if the quadrupole ordering is primary, no dipole 
moment is induced, because the doubly-degenerate eigen- 
states of Qp are the eigenstates of ±ji of Jp. 



C. RXS spectra 

Irrespective of whether the octupole or quadrupole 
ordering is realized, RXS amplitudes are generated at 
each site, according to Eq. 1)2. 3|l . They are propor- 
tional to qia2{oj){P3 -I- P4 -I- P5) for the simultaneous 
eigenstate of Tm and Qm, to qia2{uj){P3 — P4 — P5) 
for the simultaneous eigenstate of Tj^^^ and Qjn, to 
(7ia2(w)(— J-3 + P4 — P5) for the simultaneous eigenstate 
of r^^Yi ^^'^ QiTn aiid to 9iq;2(w)(-P3 - Pi + P5) for 
the simultaneous eigenstate of T^ij and Qnj. On the 
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scattering vector G = {hh£) with h + £ = odd, these 
amplitudes are summed up with a positive sign for sub- 
lattices 1 and 2 and with a negative sign for sublattices 3 
and 4. Therefore, the total RXS amplitude becomes pro- 
portional to qia2{uj)P5 for the longitudinal order, while 
they are proportional to qia2{Lo)P3 and qia2{Lo)Pi for 
the two transverse orders. Note that a similar analy- 
sis is applied to the dipole ordering. In this case, both 
the dipole and quadrupole terms contribute to the am- 
plitude. These results are summarized in Table ITTl For 
the transverse case, our present treatment could be ex- 
tended applying to the RXS spectra detected at Np M4 
edges in Uo.75Np202.'^'* In this compound, the spectra 
may be interpreted as a consequence brought about by 
the transverse type of triple-k AFO ordering driven by 
the same ordering pattern at U sites. 

Polarization dependences become particularly simple 
for G = (00£) {£ = odd) in the octupolc and quadrupole 
orderings. They are explicitly written in the scat- 
tering geometry shown in Fig. ^ as P3 = 0, P4 = 
0, P5 = (\/3/2) sin 2^^ in the a — a' channel, while 
P3 = (V3/2)cos6lcosV', P4 = (\/3/2)cos6'sinV', P5 = 
(\/3/2) sin cos 2?/; in the a — tt' channel. Figure^lshows 
the azimuthal-angle dependence of the spectra G = (003) 
in comparison with the experiment . -Qim i The experimen- 
tal data are well fitted by sin^ 2tp in the a-a' channel, 
and sin^ 9 cos^ 2iIj in the ct-tt' channel. The two trans- 
verse orders cannot reproduce the experimental curves, 
as seen from panel (b). Paixao et al. and Caciuffo et al. 
analyzed their experimental data and concluded that the 
longitudinal order gives rise to this dependencei^2i24 The 
present analysis confirms their result. Note that, based 
on a group theoretical point of view, Nikolaev and Michel 
have obtained the same resulti^ 

Now we discuss the energy profiles. In order to calcu- 
late them, we need the wavefunctions in the intermedi- 
ate state. We first evaluate the Slater integrals for the 
Coulomb interaction and the SOI parameters within the 
HFA, which are shown in Table IIIII These values are 
reduced by taking account of screening effects. The re- 
duction factors are set the same as in the ground state. 
The Hamiltonian of the intermediate state, consisting of 
the full intra-atomic Coulomb interactions between 5/- 
5/, 5/-3d and 3d-3d electrons as well as the SOI of 5/ and 
3d electrons, is represented by lOOlx (2j(i-|-l) microscopic 
states with the total angular momentum of the core hole 
jd = 3/2 and 5/2 corresponding to the M4 and the M5 
edges, respectively. Diagonalizing the Hamiltonian ma- 
trix, we obtain multiplet structures in the intermediate 
state. The 02(0;) is calculated by using Eq. (A. 8). 

The energy profile is proportional to |Q!2(aj)P in the 
octupole ordering phase. The calculated spectra around 
M4 and M5 edges are displayed with several choices of 
r values in Fig. |S1 The origin of the energy is adjusted 
such that the peak of the RXS spectrum is located at 
the experimental peak position. Since there is no reliable 
estimation for the F value, we choose three typical val- 
ues F = 0.01,0.5 and 2.0 eV. The spike-like curves with 




-45 45 90 135 

V [degree] 



FIG. 4: Azimuthal-angle dependence of the RXS spectra 
from the triple-k AFO phases with G — (003). (a) Longi- 
tudinal ordering. The open circles and filled triangles are 
experimental data, and the solid and broken lines are the cal- 
culated results, in the a-a' and a-n' channels, respectivelyi— 
(b) Transverse orderings in the ct-tt' channel. No RXS sig- 
nal is expected from these orders in the a-cr' channel. The 
solid and broken lines are the calculated results for the two 
transverse orders. 



F = 0.01 eV directly reflect the multiplet splittings of the 
intermediate states. For the M4 edge, the choice F = 0.5 
eV makes a multi-peak-structure line-shape. It merges 
into a single-peak structure around F ~ 1.0 eV. The 
choice F = 2.0 eV corresponds to one of better fittings 
with the experimental line shape. ^''■^'^ The core-level en- 
ergy is adjusted such that the calculated peak at the M4 
edge with F = 2 eV coincides with the experimental one. 
Paixao et al. reported that the line shape is well fitted 
by a Lorentzian-squared rather than a Lorentzian one^ 
As shown above, the line shape is basically determined 
by the multiplet structure, which is smeared by the life- 
time broadening. Whether it looks Lorentzian-squared 
or Lorentzian seems unimportant. As for the spectra at 
the M5 edge, their shape depends rather sensitively on 
the value of F compared to that at the M4 edge. 

The energy profile in the dipole ordering is given by 
the sum of the dipole and quadrupole terms. However, 
|Q;i(a;)p is about two orders of magnitude larger than 
|q;2(w)P. For instance, |q;i(w)P ^ 192 x |a2(w)p when 
F = 2.0 eV. Thus the dipole term usually dominates the 
quadrupole term. Although the dipole ordering is ruled 
out from experiments, we show |Q;i(a;)p in Fig. |^ as a 
reference. The peak at the M4 edge with F = 2 eV is 
at 3847.5 eV, 0.7 eV higher than the peak position of 
|q;2(w)P. Note that the spectral shape at the M5 edge 
depends on F more sensitive than that at the M4 edge. 
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TABLE II: RXS amplitudes in triple-k ordering, for G = (hhl) with h + £^ odd. 



RXS amplitude 




Longitudinal 


Transverse 1 


Transverse 2 


dipole 


-iQi(tj)(7/'o Jz Vo)(e' X e)z 
+Q2(a;){V'o|Q5|^/'o>P5(e',e) 


-iai{uj) {tPo\Jx\tPo) {e' x e)x 
+a2{uj){i'o\Q3\Tpo)P3{e',e) 


-iai{Lj){'ipo\Jy\i'o){e' x e)y 
+Q2(t^)(V'o|Q4|V'o)-P4(e',e) 


quadrupole 


a2{uj){iPo\Q5\iJo)P5{e\e) 


a2(t^)(Vo|Q3|^o>P3(e',e) 


a2 (t^) (V'o 1 Q4 1 '/'o ) -P4 (e' , e) 


octupolc 


a2{uj){iPo\Q5\i>o)P5{e',e) 


a2iu){i>o\Q3\i>o)P3{e',e) 


a2(ti')('!/'o Q4h/'o)-P4(e', e) 



Np M5 edge 



Np M4 edge 



0.00005 



0.001 





3650 



r=2.0 eV 


r=0.5 eV 




- r=0.01 eV 





r=2.0 eV 


C'Exp. G=(003) 
XExp.G=(l12) 


r-0.5 eV 




- r-0.01 eV 





3660 



3840 



3850 



(B[eV] 



tD[eV] 



FIG. 5: Energy profiles |Q:2(tt')|^ at the Np M4 (right panels) 
and Ms (left panels) edges. The lines represent the calculated 
results for F = 0.01, 0.5 and 2.0 eV, respectively, from bottom 
to top panels. In the top right panel. The open circles and 
crosses are experimental data in Np02i^^*-" The peak heights 
of them are adjusted to that of the calculated value. 



NpM, 



Np M, edge 



0.01 



0.04 



3, 





100 





3650 



r=2.o eV 
x5 


r=o.5 eV 




r=o.Ol eV 

1 







r=o.5 eV 




- r=o.Ol eV 

, ,,1 


x4 



3660 



3840 



3850 



TABLE III: Slater integrals and spin-orbit interaction pa- 
rameters in the (3d)^(5/)'' configuration within the HFA (in 
units of eV).2^ 





181.0 




29.13 




20.54 




2.158 




92.62 




2.749 


F^ 


10.39 




1.306 




59.67 


F* 


1.281 


F^ 


6.943 


G'' 


0.914 










F« 


5.017 








= 77.278 


Cs/ 


= 0.339 











Np Mj edge 


] — r=o.O eV 




— r=0.5 eV " 




r=2.0ev 


-iii^.i.i il, 





Absorption coefficients 

1 I — ^ 





3650 



Np edge 



0.5 



3660 
o)[eV] 




FIG. 7: Absorption coefficients as functions of the photon 
energy at the Np M4 (right panel) and M5 (left panel) edges. 
The bold solid and dotted lines represent the calculated re- 
sults for F = 0.5 and 2.0 eV, respectively. The vertical bars 
represent (5-functions with F = 0. 



D. Absorption coefficient 

The absorption coefRcient is proportional to 
— Im ao(ci;). We calculate q;o(w) from Eq. (A. 4) in 
the same way as in the calculation of ai{uj) and q;2(w). 
The calculated results are shown in Fig. [7| at the M4 
and M5 edges. The present calculation confirms the 
previous multiplet calculation by Lovesey et al. , in which 
— Im ao(w) has been calculated at the M4 edge for 
r = 0.7 eVm^ With increasing values of F, the multiplet 
structure merges into a single peak. The peak position 
at the M4 edge with F = 2 eV is about 0.35 eV higher 
than that in |Q;2(ct;)p. 



(o[eV] 



(n[eV] 



FIG. 6: Energy profiles |ai(w)|^ at the Np M4 (right panels) 
and M5 (left panels) edges. The lines represent the calculated 
results for F = 0.01, 0.5 and 2.0 eV, respectively, from bottom 
to top panels. 



IV. 



CONCLUDING REMARKS 



In this paper, we have studied the RXS spectra at the 
Np M4.5 edges in the triple-k multipole ordering phase 
of Np02, on the basis of a localized electron model. We 
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have derived an expression of scattering amplitudes in 
the El process, assuming that the rotational invariance 
is preserved in the intermediate states of the scattering 
process. This is a reasonable assumption when the mul- 
tiplet energy is larger than those of the CEF and the 
intersite interaction. On the basis of this expression, we 
have analyzed the RXS spectra in Np02. Assuming the 
Fg-quartet ground state, we have constructed the triple-k 
ordering ground state. The energy profiles have been cal- 
culated by taking full account of the multiplet structure 
in the intermediate state, in agreement with the experi- 
ment. 

RXS signals on multipole ordering superlattice spots 
have also been observed and analyzed at L2,3 edges 
of rare-earth metals in their compounds such as CeBg 
and DyB2C2i^"''^7r^Si?3,4n,4i.42,4.3 ^j^^ intermediate state 

is created by the transition from the 2p-core to 5c? states. 
Since the 5d states are considerably delocalized with 
forming energy bands, the assumption that the interme- 
diate state preserves the rotational invariance becomes 
less accurate. An extension of the formula is left in fu- 
ture study. 
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APPENDIX A: DERIVATION OF EQ. (1^ 

We derive a general expression of RXS amplitude un- 
der the assumption that the intermediate state keeps the 
rotational symmetry at each site. The following deriva- 
tion emphasizes the multiplet structure in the interme- 
diate state. Thereby it is more general than the previ- 
ous analyses, in which the fast collision approximation 
was adopted by replacing the multiplets with a single 
^gygl,9j44j45 ^ poxf^ of the results found in this Appendix 
were used in Ref. when we analyzed the RXS spectra 
form URu2Si2. 

Let the core hole be created at site j in the in- 
termediate state. We express the intermediate state 
as I A) = \J',M,i), where the magnitude J' and the 
magnetic quantum number M of total angular momen- 
tum (including a core-hole angular momentum) are good 
quantum numbers. To distinguish multiplets having the 
same J' value but having the different energy, we in- 
troduce the index i. Defining Maa' by Mj{e' , e,uj) = 
J2aa' <^'a'^a'Maa'{j,(^), wc rewrite Eq. (E3 as 



with 



J',M,i 

x(J',M,i|xc,'j|V'o), (Al) 



1 



huj- {Ej>^, - Eo)+iT' 



(A2) 



Assuming that the ground-state wavefunction is ex- 
pressed as a linear combination of |J, m) at each site, 

IV'o) =^cj(m) I J, m), (A3) 

and inserting this equation into Eq. (jAl|) . we obtain 
M„„,(j,c.) = c*{m)c,{m')M^f {to), , (A4) 

m,m' 

with 

J' i=l M=~J' 

X (J, TO I I J', M, i){J', M, i\xa'\J,'m'), (A5) 

where the number of the multiplets having the value 
J is denoted by Nj. We have suppressed the index 
j specifying the core-hole site. The selection rule for 
the El process confines the range of the summation 
over J' to J' = J, J ± 1. The matrix element of the 
type (J, m\xa\J', M) is analyzed by utilizing the Wigner- 
Eckart theorem for a vector operator with the use of the 
Wigner's 3j symbol;^^ 



{J\\Vi\\J') 



J' 1 J 

M fj, —TO 

(A6) 

with s±i — =F(l/\/2)(a; ± iy), so = z. The symbol 
(J||Vi|| J') denotes the reduced matrix element of the set 
of irreducible tensor operator of the first rank. Because 
of the nature of the dipole operators, M™'™ (w) ^ only 
when |to — to'I < 2. After lengthy calculation, we obtain 



M. 



a,a' (^) 



J(J+ 1) -m^ 



imai{Lj)M^^^, + ao{uj)Sa,a', (A7a) 
1 



{uj) = -/™(2TO+l)a2(cc;)(Ar"-f iM«^)„,„, 



ii/™ai(w)(Ar + iAF)a,„S (A7b) 



M. 



m+l,m 



M 



m,m+2 
a, a' 

^y-m-|-2,m 



where 



(^) 



(c.) 



-U{2m + l)a2{co)iAP 



i-fmaiiij){M- ~ iM^)„,„,, (A7c) 

a':^a2{ioKM-"'y" +iM-yUc.', (A7d) 
a':^a2{u;)iM^"~y" -iAr^Ua', (A7e) 



fn 

a'' 



x/{J -m){J + m+l), 



m+1 J 



(A8) 
(A9) 
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TABLE IV: Antisymmetric matrices, M'^ , , M'' , and sym- 
metric matrices M^"- , M^^, , M^^""''", M^'-!'^ 



Antisymmetric 



Symmetric 



0' 
1 
0-10 



NP 
0-1 



1 




1 
1 



M 



^001 

,10 

2 2 

















(;:' 


-1 


=) 


:'; 


-1 


=) 

















1 
-10 




10' 

1 




energy profiles are given by 



aoH = ^J(2J-1)F,_iH + ^J(J + 1)FjH 

+ ^(2J2 + 5J + 3)Fj+i(tj), (AlOa) 

aiH = -(2J-1)F7_iH-FjH 

+ {2J + ?.)Fj+i{uj), (AlOb) 

a2(c^) - ^[-F.7_i(c^)+Fj(c^)-F,/+i(u;)], (AlOc) 



with 



Fj,{u;) = (-)'^-^'|(J||yi||j')P^i?.(c., J'). (AH) 



and the 3 X 3 matrices, M"" , My , W , M^«, M^'^ M^^, 
A-r^'y\ and M^^'-'^' are tabulated in Table Hvl The 



Substituting Eqs. IjATp into Eq. HA4p . we obtain the final 
expression Eq. H2.3|l . 
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